Non-adiabatic spin torque investigated using thermally activated magnetic domain 

wall dynamics 
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Using transmission electron microscopy, we investigate the thermally activated motion of domain 
walls (DWs) between two positions in permalloy (Ni8oFe2o) nanowires at room temperature. We 
show that this purely thermal motion is well described by an Arrhenius law, allowing for a description 
of the DW as a quasi-particle in a ID potential landscape. By injecting small currents, the potential 
is modified, allowing for the determination of the non-adiabatic spin torque: /3t = 0.010 ± 0.004 for 
a transverse DW and /3v = 0.073 ± 0.026 for a vortex DW. The larger value is attributed to the 
higher magnetization gradients present. 
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The controlled motion of magnetic domain walls 
(DWs) by the injection of spin-polarized currents [l|, |2| 
has become an exciting field of research, since the de- 
tailed understanding of the interplay between the lo- 
cal magnetization and the spin angular momentum of 
the electron current is of scientific interest and essential 
for proposed applications [35 . Thermal effects, resulting 
from Joule heating due to the high injected current den- 
sities necessary for DW displacement, have often been 
neglected and have in general been considered as detri- 
mental in the past. However, when properly analyzed, 
thermally activated processes at temperatures even well 
below the Curie temperature can be used to provide a 
better understanding of the underlying physical mecha- 
nisms y. Current-induced DW displacement has been 
demonstrated successfully in experiments [l43|, but the 
underlying theory describing the torque exerted by spin- 
polarized conduction electrons on the local magnetiza- 
tion and, in particular, their quantitative values are still 
subject to controversial discussions [5|-l8|. 

In order to account for the effect of the current, 
two spin torque terms have been added to the Landau- 
Lifshitz or Gilbert equation the latter^o^ which reads for 
the current flowing along the x-axis 
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where Hi(t) are the effective fields, 7 = gfiB/^ is the 
gyromagnetic ratio. Si = ^u^Z/is represents the mag- 



jP^/iB/(2eMs) is an effective velocity for current den- 
sity j, P is the polarization and Mg is the saturation 
magnetization [9]. Gilbert damping is described by the 
constant a and the non-adiabaticity is characterized by 
the parameter f3. While the third term (adiabatic spin 
torque) is well understood, the physical contributions to 
the non-adiabatic torque are still the subject of scientific 
debate, and both spin- flip scattering and linear momen- 
tum transfer due to non-adiabatic transport have been 
discussed [GUs^. In particular, this torque becomes in a 
ID model an effective force [7]. While the purely adia- 
batic contribution is expected to prevail for wide DW spin 
structures, the non-adiabatic term has been proposed to 
be more effective in narrow DWs or vortices where high 
magnetization gradients occur [8|. However, the relative 
magnitude of the non-adiabatic term and its relation to 
the damping constant and to the characteristic transport 
length scale, that the DW width has to be compared to, 
are still debated and depend on the model used 
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Different methods to determine /3 have been put for- 
ward with partly contradicting results. One problem has 
been the influence of thermal effects, that stem from 
the Joule heating due to the high current densities that 
are required for DW movement and which can affect the 
physical mechanisms involved jl3|. However, a first step 
to analyze the spin torque effect using thermal effects 
was recently taken, when the thermally activated depin- 
ning processes induced by external fields, with simulta- 
neous current injection, were studied in hard magnetic 
materials [1^. DWs are pinned at artificially structured 
notches in wires, defects intrinsic to the material or de- 



fects caused by fabrication, all of which can generate at- 
tractive potential wells [2|, Il4l4l7| . Due to these strong 
pinning sites in hard magnetic materials [l4|, DWs can 
only be depinned by additional external magnetic fields 
deforming DWs prior to depinning. The resulting com- 
plicated spin structures make reliable conclusions about 
spin torque difficult without direct imaging, which has 
not been carried out so far. 

Thus, in order to use thermal effects to analyze and 
fully understand the mechanisms behind spin torque, soft 
magnetic materials with weak DW pinning are necessary 
so that thermal effects can be observed at room temper- 
ature without the application of external fields. Only 
then can very small currents sufficiently affect the mag- 
netic system without significant heating, so that the size 
of the spin torque terms [Eq. [l] can be determined using 
direct time-resolved high resolution imaging of the spin 
structure and comparison to a ID analytical description 
of thermally induced DW dynamics [71, ll8|. However, 
for DWs in soft magnetic materials such as permalloy, 
in which DWs are complicated 3D objects (see vortex 
DW spin structure [Fig. [Hd)] and transverse DW spin 
structure [Fig. [U^e)]), it is necessary to clarify if the de- 
scription based on the quasi-particle model can be used 
and if thermal depinning occurs via a single path, which 
is necessary for an analysis based on a ID model to de- 
termine the non-adiabatic spin torque. 

In this paper, we investigate in real time thermally ac- 
tivated DWs jumping between two positions in permal- 
loy (Ni8oFe2o) nanowires using transmission electron mi- 
croscopy. This motion occurs at room temperature and 
is of purely thermal origin, in the absence of an external 
magnetic field or injected electron current. The distri- 
bution of the dwell times for which a DW stays at each 
of the two positions is well described by an Arrhenius 
law with the DW described as a quasi-particle in a ID 
potential, with two metastable states separated by an 
energy barrier. By modifying the local potential using 
a small constant current with no significant heating and 
by analyzing the dwell times, we infer values for j3 for a 
transverse and a vortex DW. We find that there are sig- 
nificant differences between the values determined for /3 
for different DW types, which can be attributed to their 
different spin structures. 

Fig. [B^a) shows a transmission electron microscopy im- 
age of a permalloy zigzag wire on a silicon nitride mem- 
brane. These structures (wire width: 150-500 nm, wire 
thickness: 8-20 nm) were fabricated by electron-beam 
lithography and a lift-off procedure [19|. At both ends, 
the line structures are contacted by Au pads, allowing 
for the injection of currents. To improve heat dissipa- 
tion, the membrane substrates can be additionally back 
coated with 30 nm of Aluminum. 

Transmission electron microscopy is a particularly suit- 
able technique to observe and understand the full ther- 
mally activated stochastic dynamics of DWs as it pro- 
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FIG. 1: (a) Bright field TEM image of a permalloy zigzag 
line structure. The arrows indicate the directions of magne- 
tization and the positions of nucleated DWs are marked by 
circles, (b) and (c) Fresnel images of thermally activated vor- 
tex and transverse DW jumping between two pinning sites, 
(d) and (e) Electron holography magnetic induction maps of 
vortex and transverse DW in the same wires showing detailed 
spin structures indicated by the color code and the arrows. 



vides high spatial (< 10 nm) and time resolution (< 0.05 
s) for single shot measurements. Here, the Fresnel mode 
of Lorentz microscopy is used to identify the positions of 
the DWs from the dark or bright contrast that results 
from the deflection of electrons by the Lorentz force [20| . 
For the investigation of detailed DW spin structures, we 
use off- axis electron holography [21]. 

The sample is initially magnetized using a field of '^ 1 
T in a direction perpendicular to the wires. Tilting the 
sample by 30° results in the nucleation of DWs in the 
kinks of the wires after relaxing the field (see schematic 
in Fig. [B^a)). We find either vortex [Fig. [TJb) and (d)] 
or transverse [Fig. [TJc) and (e)] DWs depending on the 
wire geometry [15]. 

We first analyze a transverse DW (sample 1: wire 
width 495 nm, wire thickness 10 nm) jumping between 
two pinning sites (separation 318±32 nm) without the 
influence of an external magnetic field or current. This 
pure thermal DW motion occurs at room temperature 
since no significant heating effects due to the electron 
beam of the transmission electron microscope were ob- 
served (The sample temperature is monitored by mea- 
suring the resistance [22|). After acquiring 7500 Lorentz 
images of the DW (DW width: 590 ±10 nm) in the movie 
mode, we identify the DW positions and investigate the 
dwell times ri and T2 that the DW stays at each of the 
two positions. 

The exponential decay of the dwell time distributions 
[Fig. [2], and in particular the fact that it is well approx- 
imated by a single exponential function [23|, means that 
the DW can be considered as a quasi-particle moving in 
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FIG. 2: Normalized distribution of dwell times for both pin- 
ning sites (n and T2) of a thermally activated transverse DW 
in a permalloy wire (sample 1) fitted with a single exponential 
function. 



a ID potential between two metastable states separated 
by an energy barrier, with a single transition path across 
it (shown schematically in Fig.[3fa)) [14]. In general, the 
pinning strength is not equal at the two positions, yield- 
ing two depths for the potential wells with two different 
energy barriers eo,i and 60,2- The probability of a tran- 
sition to the neighboring state is only determined by the 
ratio between the corresponding energy barrier and the 
thermal energy that excites the transition. This depen- 
dence results in finite dwell times t\ and T2 , for which the 
DW stays at each of the two local metastable states be- 
fore it passes to the neighboring minimum of the pinning 
potential. 

We carry out the same investigation for a vortex DW 
(sample 2: wire width 395 nm, wire thickness 17 nm). 
The two pinning positions (separation 20 ± 4 nm) of the 
DW (DW width 363 ±10 nm) are again extracted from 
Lorentz [Fig.[T]^b)] and holography images [Fig.[T]^d)] and 
we find that the complete DW is displaced. The dwell 
times show the same exponential decay as observed above 
for the transverse DW suggesting that the DW transi- 
tions can again be described by an Arrhenius law. 

Since the dwell times of both DW types obey an Ar- 
rhenius law, we use in the following a ID quasi-particle 
description [7| to analyze the influence of a constant di- 
rect current on the dynamics of the thermally activated 
DWs. In this model, the non-adiabatic torque acts as a 
force on the DW, Fdw = —N?il3ux/\^ where N is the 
number of spins in the DW, iix is an effective spin cur- 
rent drift velocity and A is the DW width [7, 24, 25]. This 
force affects the ID potential [Fig. [3l^b)] and the current 
dependent shift in the energy Cj is: 
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FIG. 3: (a) Double well potential for a ID description 
of thermally activated DWs jumping between two positions 
(schematically shown for a transverse DW) . For the case with- 
out a current, the energy barriers consist of current inde- 
pendent contributions eo,i and eo,2- The difference between 
the two barriers is Aeo, the distance between the two pin- 
ning sites is xq. (b) An additional force on the DW results 
from the injection of a spin-polarized current due to the non- 
adiabatic spin torque. The energy barriers now consist of 
current independent eo,i and current dependent contributions 
ej^i (i = 1, 2). (c) An example of measurements of the time- 
resolved thermal movement. 



sectional area of the DW, a^ the volume of the unit cell 
and the spin polarization (P = 0.37 [26|). The charac- 
teristic dwell times, for which the DW stays at each of 
the two positions are described by the Arrhenius law [j] : 



II = !M.e(^o>i+^i>i-(^o>2+e„2))/(feBT) 
T2 Too 



(3) 



Since the small current densities only induce a linear 
contribution to the current independent pinning poten- 
tial (Eq. [2j), the curvature of the potential at the two 
metastable states, which is given by the second deriva- 
tive, is not significantly affected by the applied currents. 
Thus, the attempt frequencies 1/tq,i and 1/to,2 deter- 
mined by the potential curvature [45 are independent of 
the current. From Eq. [2]and Eq. [3l we obtain: 
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By measuring dwell times as a function of current density 
and fitting these with Eq. HI we obtain values for f3. 

For a transverse DW (sample 1), we measure the de- 
pendence of dwell times on the intensities and directions 
of constant direct currents. For each current density, 
7500 Lorentz images are acquired in order to obtain suf- 
ficient statistics when extracting the positions of the DW 
(see example in Fig. [3l^c)). The resultant dwell times for 
both metastable states [Fig. ID^a)] are fitted to the ID 
model (Eq. |4]) in order to calculate /3 from the slope of 
the fit. From the Lorentz images and the electron holo- 
grams, values for the jump distance xq and the DW width 
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FIG. 4: Measurements of dwell times n and T2 plotted as 
a function of applied direct current for (a) transverse and a 
(b) vortex DWs. The non-adiabatic coefficient (5 is calculated 
from the slope of the fit with the ID model (Eq. [3|). 



values is attributed to the high magnetization gradients 
at the vortex core. The comparatively large value for 
/3v means that high DW velocities (> 100 m/s) should 
be possible for vortex walls, which bodes well for fast 
switching devices based on DW motion. 
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A are derived and from the fit in Fig. HJ^a) we obtain 
Pt = 0.010 ± 0.004 as result. This value is similar to the 
value of the Gilbert dammag constant measured for our 
permalloy, apy ^ 0.008 [27|. 

For the vortex DW (sample 2), positions and dwell 
times are again extracted from Lorentz images, as shown 
in Fig. |li;b), to obtain a value for p^ = 0.073 ± 0.026, 
which is considerably larger than the Gilbert damping 
constant. This measurement results in a larger non- 
adiabatic coefficient for the vortex DW than for the 
transverse DW (/3t ^ 1.3 • apy for the transverse DW; 
/3v ^ 9.2 • apy for the vortex DW). 

The fact that /3v ^ A can be explained by using 
a more sophisticated description of the non-adiabatic 
torque, which includes two contributions /3v = /3" + /3^ 
with /3" accounting for non-adiabatic transport and j3^ 
for spin relaxation due to spin- flip scattering |28|. Since 
the magnetization inside a transverse DW varies slowly, 
the corresponding magnetization gradients are small, and 
non-adiabatic contributions to the spin torque effect are 
mainly the result of spin- flip scattering due to impurities 
and spin-orbit coupling [7]. Therefore, j3t ^ P^ is ex- 
pected for a transverse DW in permalloy and values of the 
order of the damping constant a are predicted 12|, |29| . 
For a vortex DW in the investigated permalloy nanowires, 
the magnetization turns out-of-plane at the vortex core 
and a large magnetization gradient occurs. Thus, non- 
adiabatic transport due to reflected conduction electrons 
becomes more significant and /3" is larger [8[, while the 
contribution due to spin relaxation j3^ is still present. As- 
suming that the contribution from spin relaxation j3^ is 
similar for transverse and vortex DWs as it is intrinsic to 
the material, we can derive the pure contribution of non- 
adiabatic transport /3" = /3v - A = 0.063 ± 0.030. We 
note that these results are deduced from analysis based 
on a ID model, which nevertheless seems to be a reason- 
able assumption given the good fit to the Arrhenius law 
in Fig. El 

In conclusion, we derive p from thermally activated 
vortex and a transverse DW obtaining /3v = 0.073 ± 
0.026 > A = 0.010±0.004. The difference between these 
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